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Abstract. We consider different generalizations of tlie Euler formula and dis- 
cuss the properties of the associated trigonometric functions. The problem is 
analyzed from different points of view and it is shown that it can be formulated 
in a natural way in algebraic and geometric terms. 



1. Introduction 

The circular trigonometric functions have long been recognized as the charac- 
teristics functions associated with the exponentiation of imaginary numbers and 
more in general of matrices. Within such a framework the notion of trigonometry 
is not unique and we can consider different "trigonometries" whose common thread 
can be traced back to some underlying algebraic structure. This point of view has 
inspired the researches of Ref. where a systematic group theoretic derivation 
of different formulations of trigonometries has been undertaken. The need for a 
wider conception of trigonometry has been pointed out in Refs. [H [3], where the 
trigonometry in Minkowskian spacetime has been studied and exposed in depth. 

Circular functions and complex numbers are intimately connected and the most 
representative link is provided by the Euler identity 



(1.1) 



id 



The properties of the circular functions can be derived, using a purely algebraic 
point of view, from eq. (jl.ip and from the cyclical properties of the circular imagi- 
nary unit I, satisfying the identity — —1. Furthermore, the differential equations 
satisfied by the ordinary trigonometric functions are just a consequence of the fact 
that the exponential function is an eigenfunction of the derivative operator, while 
the addition formulae are a consequence of the semi-group property of the expo- 
nential. 

The introduction of a hyperbolic "imaginary" unit such that = 1 allows 
the definition of the hyperbolic functions, which satisfy the identity 

(1.2) e^^ = cosh6l-^jsinh6l. 

The theory of the trigonometric functions (in their wider form, circular and 
hyperbolic) can be developed independently of any geometrical interpretation and, 
within such a context, the distinction between the units i and j does not appear 
particular meaningful: both can be viewed as solution of a second degree algebraic 
equation. We can therefore consider the following realization of a more general 
"imaginary" unit [5] 

(1.3) h^^a + bh (a,6eR). 
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By interpreting this equation as a second degree algebraic equation, its solutions 

(1.4) ^± = ^^ (A = fe2 + 4a) 
can be used to introduce the following identities 

(1.5) e''^'^ = c{e) + h±s{e), 

that can be viewed as a generalization of the Euler formula. We have indicated 
with C and S two functions of the variable 6 that play the role of the trigonometric 
functions in eqs. (jl.ip and ()1.2|) . For this reason we refer to them as trigonometric- 
like functions (TLF) . It is worth stressing that in writing eq. ()1.5p we have assumed 
that h± are conjugated imaginary units, i.e., h'^ — h^, as in the case of circular 
and hyperbolic units. 

By combining the previous equations, we can identify the functions C and S as 

(1-^^ [s{e) j"VA I 1 -1 j U'^-' 

The differential equations satisfied by this family of functions is obtained by keeping 
the derivative with respect to 9 of both sides of eq. (|1.5|) with respect to 9 and, by 
using eq. (|1.3p . we end up with (we use h without any index because the identity 
depends only on the property (II. 3p of the new imaginary unit) 

C'{9) + hS'i9)^hC{9) + {a + bh)S{9), 

that, after equating the coefficients of terms in h, yields the following first-order 
differential equation 

(1.7) C'{9)^aS{9) S'{9) ^C{9)+bS{9). 

Moreover, it is easy to obtain the addition formulae [6] 



(1.8) 



Ci2 A / Ci aSi \ { C2 



where we put 

fk = f{9k) fki = f{9k + 9i) if = C,S;k,l ^1,2). 



It is worth noting that the functions C and S derived from eq. (11.51) can be 
expressed in terms of ordinary trigonometric functions (circular or hyperbolic), 
and therefore they cannot be therefore considered a new set of functions. How- 
ever, if considered independently from the existence of their classical counterparts, 
they provide a tool having global properties (e.g., the addition formulae and the 
differential equations they satisfy) which make them very useful in applications, 
as, for example, the evolution of multilevel quantum systems (as discussed in the 
concluding section) , or to treat problems of analytical mechanics [3 [7] . 

The following general result holds: Given a, /3 G K, it is always possible build a 
set of cosine- and sine-like functions according to the following prescription 

(1-9) C{9) = ^ ^ 3(9) = , 

a — p a — p 

whose proof is straightforward and is achieved by treating a,f3 as the independent 
roots of a second degree equation. 
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2. Trigonometries and 2x2 matrices 

The notion of imaginary units exploited in the previous section can be made 
less abstract by noting that both circular and hyperbolic imaginary units can be 
expressed in terms of matrices, namelj0 

in such a way that 

(2.2) = f = i , 

with 1 the unit 2x2 matrix. As well known, the exponentiation of t (the symplectic 
matrix) can be viewed as the generator of rotations in a plane, while the exponen- 
tiation of J can be associated with non Euclidean rotations, and, more in general, 
with the Lorentz transformations viewed as rotations in the Minkowski plane. 

By introducing the vector z{d) ~ {C (9) , S (9))'^ , eq. (jl.7|) can be rewritten as 
follows 

(2.3) ^z{9) = hzi9), ^ = ( ? 6 

with the matrix h satisfying the identity (jl.3|) . According to this equation the TLF 
are associated with the non-unitary "evolution operator" 

(2.4) U{9)^e^^, 

The use of standard algebraic methods |8J allows to derive the following expressions 
for the entries of the matrix associated to this operator 

( V. 

(2.5) U22i9) = 

Ui2i9) ^. 2 . / VA 

= U2i{9) = -j= smh — 

a VA Y 2 

It is worth noting that the functions C and S can be viewed as independent linear 
combinations of the entries of the matrix (|2.4|) . or, what is the same, that, according 
to eq. ()1.8|) . the evolution matrix can be written as 




(2.6) uie)=(^^'^ «^(^) 



s{e) c{e) + bs{e) 

The previous results can be extended to any generic 2x2 matrix 

In fact, the Cayley-Hamilton theorem [19; allows to prove the identity 
(2.8) = -(detrn) i + (trm)m 



^In the following, a superimposed "hat" denotes a matrix. 
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that is a generahzation of eq. (|1.3p . and, therefore, ah the functions that satisfy 
the differential equation 

(2.9) ^z{e)^mz{e), 

are TLF. The use of eq. (|1.3I) with rh in place of h, allows to write the entries of 
the evolution matrix associated to eq. (12. 9p as follows 

- 1, c5(0) C((?)+d5(0) ) ' 

that, compared with the expression of the matrix entries algebrically obtained 

(A = (tr rhf - 4 dot to = (a - d)^ + 4 6 c) , 

allows to obtain the explicit form for these TLF. As a consequence of eq. (|2.10p . 
the following identity holds 

(2.11) exp(6ltrTO) = C^{e) + (tr m) C(6l) 5(6*) + (det to) 5^(61) , 

whose relevant geometrical and physical meaning will be discussed in the concluding 
section. 

3. Higher order trigonometries 

According to the analysis developed so far, what we have called trigonometry 
is the by-product of the formalism associated with the exponentiation of 2 x 2 
matrices. Higher order trigonometries can accordingly be considered along with 
the exponentiation of matrices with larger dimensionalities. 

By following Ref. ^Zi we propose a further generalization of the Euler identity, 
by setting 

(3.1) c^' =Aa{e)+i^A^{e)+7fA2{e), rf = a^ + T^a^+rf a2. 

The expression of the functions Aj can be found adopting exactly the same proce- 
dure as before, albeit slightly algebraically more cumbersome. By setting 

A{e) = {A,{9),A,{9),A2{9)f , Y{9) = (e"^ ^ e"^ ^ e''^ , 

with rjj (j — 1,2,3) roots of the cubic equation in eq. p.ip . we obtain 

(3.2) A{9) = V-^Y{9), 



We adopt the symbols Aj(j = 0, 1, 2) for the trigonometric- Hke functions because more con- 
venient and easily amenable for the n-dimensional extensions. 
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where V is the Vandermonde matrix 

/ 1 VI 

V = \ 1 V2 

V 1 '73 

Since a triple of numbers can always be associated with a cubic equation, it is clear 
that we can construct a set of third order TLF for each of these triple following the 
prescription (|3.2|) . 

The differential equations satisfied by the third order TLF are derived by means 
of the same procedure as before, which yields 

(3.3) ^Aie)^fiAi9) 
with 

/ ao 

(3.4) 77= 1 ai 

V 1 a2 

The evolution matrix associated to the previous Cauchy problem can be written as 

(3.5) U{0) ^e'^ ^Ao{9)+f,A,{e)+ ff A^iO) , 
and, in terms of the functions Aj, reads 

(Aq ao A2 ao Ai + ao 02 A2 

Ai Aq + ai A2 ai Ai + (ao + ai 02) A2 
A2 Ai+a2A2 Ao + a2Ai + {ai+al)A2 

This equation can be used to compute the determinant of the evolution matrix in 
terms of the two-variable TLF, and from the relation 

(3.7) detC/ = exp(6'tr77) 
we obtain 

e°'* = Al+aoAl+ai A^ - (300 + 0102)^0^1^2 

(3.8) + (2 01 + a^) A2 Al + (o^ + 2 oo 02) Aq Aj 

— oi Af Aq + Oo 02 Af Aq - Oo oi Ai Al 

that should be viewed as the fundamental identity for the third order trigonometry. 

As a particular case of third order TLF, let us consider the case of the Eisenstein 
numbers [10]. They satisfy the properties 

(3.9) uj^ = l, 

i.e., can be obtained letting oo = 0, oi = 02 = — 1 and 77"^ = 1 in the second of eq. 
(|3.ip . Following the treatment of sec. ([1]), the associated C and 5* functions are 
easily obtained 

■73 



0(9) = e 



- p-e/2 



cos 





(3.10) 

It is interesting to link the above functions to the so called pseudo-hyperbolic func- 
tions (PHF), discussed on the eve of the seventies of the last century [TT], and more 
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recently discussed within the framework of an algebraic formalism in Refs. [12j and 
[l3] . According to eqs. p.9p we find that 



(3.11) 
with 



g3 n+fc 



^ 3n 

n=0 ^ 



A:)! 



(^^ = 0,1,2), 



Furthermore, from the condition 

(3.12) e"'' = C(6i)+wS'(6i), 
by using eqs. (13. 9p and (|3.11l) . we obtain 

(3.13) C{B) = eo{e)~e2{e), 3(6) ^ ei{0) - e^iO) . 

The extension of these results to the rt-dimensional case is straightforward and 
computationally simplified by the remark that the structure of the matrices (12. 3p 
and p.4[) is that of a companion matrix We remind that, given a n x n matrix 
Q with characteristic polynomials 



(3.14) 



k=0 



its companion matrix is defined as 



(Po - 1) , 
\ 

/ 

and has the same eigenvalues of Q. Introduced the n-dimensional vector e 
(1, 0, • • • , O)"'", it is easy to show [9 that the following identity holds 







( 


• 


• 








1 


• 


• 


~Pn-l 


(3.15) 


L = 





1 • 


• 


-Pn-2 






I 


• 


• 1 


~Pl 



(3.16) 



e = 



ji-i 

E 

fc=0 



This result allows to conclude that 



(3.17) 



and, since e^^ can always be written in terms of n suitable TLF, wc can state that 
the exponentiation of any diagonalizable n x n matrix can be written in term of 
these functions. 

Before concluding this section, let us come back to eq. p.5|) and consider its 
extension to the case of two variables 



(3.18) 



[/(0,0) 



+ fiAiid,(f>) + f A2i 



Such an extension implies that the functions Ak are actually two variable functions. 
Their forms are obtained quite straightforwardly by noting that the "evolution" 
matrix (|3.18p and r) are diagonalized by the same transformation. The functions 



are therefore obtained by performing the substitution e 



Vj I 



, r,j e+7j^ 4, 
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in eq. p.2p . Let us note that in this case, the determinant of the matrix [7(6', (/>) is 
given by 

(3.19) detU{e,(j)) :^exp{eti-fi + (f>tyif) 
where 

3 3 

(3.20) tr ?7 = ^ rjj = a2 , tr 77^ = ^ 77^ = 03 + 2 ai . 

and therefore the only modification in eq. (j3.8|) consists in the multiphcation of its 
Ihs by the term e(''2+2ai)0. 

4. Concluding remarks 

As the usual trigonometric functions (circular or hyperbolic) , also the TLFs have 
a geometric interpretation. Let us consider the identity given in eq. (j2.11l) . In terms 
of the variables 

X = exp ^— - tr rh^ C{0) , y = exp 2 ^ ^^-^ ' 

it can be rewritten as follows 

(4.1) + {tr m) X y + {det m) ^ 1 . 
The determinant of this quadratic form is 

r (trTO)2 , , ^ , (a - d)^ 

(4.2) ^ ^' - (det m) = — -^-!-+bc 

and its sign determine the geometric interpretation. For 5 > Q the points with 
coordinates x, y, lie on an arc of hyperbola, and reduce to the ordinary hyperbolic 
functions when referred to the axes rotated with respect to x, y by an angle 

(4.3) X = ^ tan^i 



2 \ad-bc- 1^ 

In the case ^ < the quadratic form is represented by an ellipse and the same 
axis rotation reduce these functions to the circular ones. Finally, for 5 = the x, y 
functions does not make too much sense because the conic is degenerate. 

As an example of third order TLF, we consider the particularly important case 
in which the matrix 77 has the following form 

-V3 1^2 

(4.4) r}= I -vi 

-V2 Vl 



This matrix is not of the form p.4p and therefore the associated functions Ak 
cannot be considered canonical. By using eq. (13.21) it is easy to show thatH 



(4.5) e'^^ = i + eiiinc{ve)ff+-e'^fimc^{-6)ff i^^^Yl 

\ 1=1 J 

A formalism based on the above concepts has been exploited in Ref. [15] as a 
tool to work with the matrices describing the dynamics of the quarks and neutrinos 



^sincx = sinx/x. 
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mixing. In this context, the matrix r} is recognized as the generatrix of the Cabibbo- 
Kobayashi-Maskawa [16] and Pontecorvo-Maki-Nagakawa-Sakata [17] matricefl 

In sec. [T] we have mentioned the addition theorems, which are derived as a 
consequence of the semi-group property of the exponential function. But what 
about e^(''i+''2\ with hi and ft,2 both satisfying eq. (|1.3p ? Under this hypothesis, 
taking into account that the "numbers" /11.2 in general do not commute, one has 

(4.6) {hi + h2f = a + bihi+b2h2 + hih2 + /i2 hi , 
that, unless 

(4.7) hi h2 + h2hi= hi hi + 62 /ii , 

cannot be put in the form (jl.3p . In the 2-dimensional case the problem can be 
overcome by noting that we can always write 

3 

Ci Gi 



(4.8) hi + h2 = l + ^i 



i=l 



where CTi are the Pauli matrices, and therefore we get 

3 

(4.9) {hi + h2f ^a + 2{hi + h2) a=-l + ^c^ 

1=1 

i.e., the same structure of eq. (|1.3p . This implies that, analogously to eq. (|1.5p . we 
can write 

(4.10) Qe{hr+h2) ^ ^(g,) _^ ^.(g,) 

For higher order representation we can use different realizations of the Clifford 
numbers, e.g. the Dirac matrices in the 4-dimensional case. 

The link of the above procedure to the Hermite-heat polynomials and integral 
transforms of the Airy type will be discussed in a forthcoming investigation, where 
we will also analyze the properties of "complex" functions for which the following 
decomposition holds 

(4.11) f{x + hy)^u{x,y) + hv{x,y) 

where the functions u and v satisfy the Cauchy-Riemann condition 

(4.12) dyw{x,y) ^ hd.j;w{x,y) (w=(m,w)^) 

that, taking into account eq. (|2.3p . reflects in the following partial differential 
equation to be satisfied by the functions u and v 

(4.13) {dl~adl-bdly)v.^O. 

Therefore, the functions u and v can be considered a generalization of the ordinary 
harmonic functions. 

Furthermore, if the function / admits a Fourier transform, we can write 

f{x + hy) - -= / dfce''=(-+''^)/(fc) 

V2 TT J-00 

1 p°° 

(4.14) ^ -= Ake'^'-" f{k)[C{iky) + hS{iky)] 



/2 7r 



■^These matrices have complex entries, while all the elements of matrix 1 14.411 are real. 
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thus getting for the functions u and v the foUowing integral representation 

The previous relations can be viewed as an extension of the concept of the ordinary 
Fourier transform. A further element which can be brought to the discussion is the 
possibility of extending the realm of special functions by introducing new families 
of Bessel functions. We can indeed use the following extension of the Jacobi- Anger 
generating function 

oo 

(4.16) e^^W= e("'"+''")^i3„,,„(a;) 

m,n— — cxD 

to introduce the two-index Bessel functions i?„i^„(x), with S{9) given by eq. (|1.9|) . 
The recurrence relations satisfied by B„iji{x) can easily be obtained from their 
definition itself. By keeping the derivative of both sides of eq. (14.16^ with respect 
to X, and equating the like-power coefficients in e"™^ and e'^"^, we obtain 

(4.17) (a - /3) Bm.n{x) = Bm-l,n{x) - B„i^n-l{x) , 

da; 

while by deriving with respect to 9 we find 

(4.18) {am + l3n)Bm,nix) = — ^ [a B^^i,n{x) - l3 Bjn.n~i{x)] . 

a — p 

It is evident that an analogous function can be associated with the cos-like part- 
ner of the S{9) function, and that the same procedure can be exploited for the 
higher order trigonometric functions. This is a fairly promising way to improve the 
theory of Bessel like functions. These type of problems, along with those relevant 
to properties of the complex functions (|4.11|) will be considered in a forthcoming 
investigation, where we will study more carefully the problem associated with the 
non-commutativity and the exponential disentanglement. 
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